Maple 2018 Questions and Posts

Computation of left and right eigenvectors...

Asked by: Elisha 50

June 17 2025

0 0

Someone should please help me compute the left and right eigenvectors of the system below. The purpose is to compute values for 'a' and 'b' in the bifurcation formula.

Thank you

>

(1)

>

(2)

>

(3)

>

diff(S(t), t) := `Λ__p`-(`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p]+µ__C)*S+`ω__B`*I__B

Lambda__p-(`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p]+Âµ__C)*S+omega__B*I__B

(4)

>

diff(I__B(t), t) := `#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S/N[p]-`ω__B`*I__B-(`σ__B`+µ__C)*I__B

`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("θ",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S/N[p]-omega__B*I__B-(sigma__B+Âµ__C)*I__B

(5)

>

(6)

>

diff(S__A(t), t) := `Λ__A`-(µ__A+`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p])*S__A+`δ__A`*I__A

Lambda__A-(Âµ__A+`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`/N[p])*S__A+delta__A*I__A

(7)

>

diff(I__A(t), t) := `#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S__A/N[p]-(µ__A+`δ__A`)*I__A

`#mrow(mi("ϕ",fontstyle = "normal"),mo("⋅"),msub(mi("α",fontstyle = "normal"),mi("B")),mo("⋅"),msub(mi("I"),mi("B")))`*S__A/N[p]-(Âµ__A+delta__A)*I__A

(8)

>

(9)

Download CBD2.mw

Simplifying a nested radical...

Asked by: Kitonum 21430

June 05 2025

3 2

It is not difficult to manually check the validity of the identity sqrt(x+2*sqrt(x-1)) = sqrt(x-1)+1 , which is true at least for all x>=1 . I don't know of any Maple command that directly makes this simplification:

expr:=sqrt(x+2*sqrt(x-1));
simplify(expr);
simplify(expr) assuming x>=1;

It was possible to simplify it only by making a variable substitution x=y+1 (and then a reverse substitution y=x-1 ) :

expr1:=subs(x=y+1,expr);
simplify(expr1);
subs(y=x-1, %);

By the way, the CAS Mathematica also cannot cope with simplifying the expression expr .

Pde Question How to impliment the this equations i...

Asked by: KIRAN SAJJAN 40

May 16 2025

0 11

How to solve these pde equations in maple to get the similar type graphs.

Ode equations we can solve directly but these equations are pde .

in the article they have solved the governing equations by series solution?

can we solve these equations in maple also by series solution or any other method is there to solve these equations

Issue in Maple 2018.2...

Asked by: Kitonum 21430

December 17 2024

1 4

is(abs(x)=max(x,-x)) assuming real;

# FAIL

I wonder if this will work in newer versions of Maple?

How can I simplify this relation?...

Asked by: mary120 75

December 13 2024

0 3

Hi,
How can I simplify this relation(See uploaded .mw file)?
For example, the second term is simplified as:

deltae*(1-phi0/(kappa-3/2))^(-kappa+1/2)+(1/2)*deltab*(1-sqrt(2)*sqrt(1/(m*ub^2))*sqrt(-phi0));

di1.mw

conflict in equation names in Section and main fil...

Asked by: PinheiroRJ 5

November 14 2024

1 6

How to insert a file inside another file, creating a session, without equations with the same name conflicting?

How to solve this equation?...

Asked by: mary120 75

November 03 2024

0 3

Hi,
I have an equation and I want to solve it parametrically to find x , but I couldn't do that with "solve" command. (I know x should be real and positive). What should I do?
Root_of.mw

Apply Keller Box method for ODE...

Asked by: KIRAN SAJJAN 40

February 20 2024

0 2

Dear sir ,

I have implemented Dsolve method the code was executed, but i need to apply Kellor Box method to solve the ODES

Please can any one help how to implement?

because there is no post regarding the Kellor box method.

>

for j to numelems(KpVals) do Ans[j] := dsolve(eval([OdeSys, Cond], Kp = KpVals[j]), numeric, output = listprocedure) end do

>

with(plots):
cols := [red, blue, black,green]:

plotA:= display
  ( [ seq
      ( odeplot
        ( Ans[k],[eta,(f(eta))],
          eta=0..1,
          color=cols[k]
        ),
        k=1..numelems(KpVals)
      )
    ],linestyle = "solid",
    'axes'= 'boxed',labels=[eta,'f(eta)'],labelfont=[TIMES,BOLD,16]
  );

>

with(plots):
  cols := [red, blue, black,green]:

plotB:= display( [ seq( odeplot
        ( Ans[k],[eta,Theta(eta)],
          eta=0..1,
          color=cols[k]
        ),
        k=1..numelems(KpVals)
      )
    ],linestyle = "solid",
    'axes'= 'boxed',labels=[eta,'Phi(eta)'],labelfont=[TIMES,BOLD,16]
  );

>

with(plots):
  cols := [red, blue, black,green]:

plotC:= display( [ seq( odeplot
        ( Ans[k],[eta,Phi(eta)],
          eta=0..1,
          color=cols[k]
        ),
        k=1..numelems(KpVals)
      )
    ],linestyle = "solid",
    'axes'= 'boxed',labels=[eta,'Phi(eta)'],labelfont=[TIMES,BOLD,16]
  );

>

with(plots):
cols := [red, blue, black,green]:

plotA:= display
  ( [ seq
      ( odeplot
        ( Ans[k],[eta,(diff(f(eta),eta))],
          eta=0..1,
          color=cols[k]
        ),
        k=1..numelems(KpVals)
      )
    ],linestyle = "solid",
    'axes'= 'boxed',labels=[eta,"f '(eta)"],labelfont=[TIMES,BOLD,16]
  );

>

Download kellor_box_method.mw

How to solve this Eq. by Maple?...

Asked by: raj2018 55

February 14 2024

0 2

Hi,
I have a problem and I haven't been able to solve it yet. I want to solve an ordinary diffrential equation similar to
(dphi/dxi)^2+2*V(phi)=0
and plot phi versus xi for a the following conditions:
1) V(phi)=dphi/dxi=0 at (phi=0,phi_m) and
2) dV(phi)/dphi=0 at phi=phi_m and
3) d^2V(phi)/dphi^2=0 at both phi=0 and phi=phi_m.
How can I do this by Maple?(see the attached file)
w1.mw

Ode equations solution ...

Asked by: SHIVAS 35

December 15 2023

0 4

How to solve this type of ode in maple

need the value of S, Q, E

unable to plot 3d and 2d plots...

Asked by: KIRAN SAJJAN 40

October 05 2023

0 3

I am unable to draw both 3d plots sowing error please help me to solve

>

restart:

>

B1 := 1+2.5*Px+6.2*Px^2; B2 := 1+13.5*Px+904.4*Px^2; B3 := 1+37.1*Px+612.6*Px^2; B4 := (ks1+2*kf-2*Px*(kf-ks1))/(ks1+2*kf+Px*(kf-ks1)); B5 := (ks2+3.9*kf-3.9*Px*(kf-ks2))/(ks2+3.9*kf+Px*(kf-ks2)); B6 := (ks3+4.7*kf-4.7*Px*(kf-ks3))/(ks3+4.7*kf+Px*(kf-ks3))

>

a5 := 1+3*((p1*sigma1+p2*sigma2+p3*sigma3)/sigmaf-p1-p2-p3)/(2+(p1*sigma1+p2*sigma2+p3*sigma3)/((p1+p2+p3)*sigmaf)-((p1*sigma1+p2*sigma2+p3*sigma3)/sigmaf-p1-p2-p3))

>

ODE:=[(a2+K)*(diff(U0(eta), eta, eta))/a1-Ra*(diff(U0(eta), eta))+lambda0/a1-a5*M1^2*U0(eta)/a1+K*(diff(N0(eta), eta))/a1+la*Ra*Theta0(eta)*(1+Qc*Theta0(eta)), (a2+K)*(diff(U1(eta), eta, eta))/a1-H^2*l1*U1(eta)-Ra*(diff(U1(eta), eta))+lambda1/a1-a5*M1^2*U1(eta)/a1+K*(diff(N1(eta), eta))/a1+la*Ra*(Theta1(eta))(1+2*Qc*Theta0(eta)), diff(N0(eta), eta, eta)-Ra*a1*Pj*(diff(N0(eta), eta))-2*n1*N0(eta)-n1*(diff(U0(eta), eta)), diff(N1(eta), eta, eta)-Ra*a1*Pj*(diff(N1(eta), eta))-2*n1*N1(eta)-n1*(diff(U1(eta), eta))-H^2*a1*Pj*l1*N1(eta), (a4/(a3*Pr)-delta*Ra^2/H^2+4*Rd*(1+(Tp-1)^3*Theta0(eta)^3+3*(Tp-1)^2*Theta0(eta)^2+(3*(Tp-1))*Theta0(eta))/(3*a3*Pr))*(diff(Theta0(eta), eta, eta))-Ra*(diff(Theta0(eta), eta))+a5*Ec*M1^2*U0(eta)^2/a3+(a2+K)*Ec*(diff(U0(eta), eta))^2/a1+Q*Theta0(eta)/a3+4*(diff(Theta0(eta), eta))^2*Rd*(3*(Tp-1)+6*(Tp-1)^2*Theta0(eta)+3*(Tp-1)^3*Theta0(eta)^2)/(3*a3*Pr), (a4/(a3*Pr)-delta*Ra^2/H^2+4*Rd*(1+(Tp-1)^3*Theta0(eta)^3+3*(Tp-1)^2*Theta0(eta)^2+(3*(Tp-1))*Theta0(eta))/(3*a3*Pr))*(diff(Theta1(eta), eta, eta))-(H^2*l1+2*Ra*delta*l1+Ra)*(diff(Theta1(eta), eta))+(Q/a3-delta*H^2*l1^2)*Theta1(eta)+2*(a2+K)*Ec*(diff(U0(eta), eta))*(diff(U1(eta), eta))/a1+2*a5*Ec*M^2*U0(eta)*U1(eta)/a3+4*(diff(Theta0(eta), eta, eta))*Theta1(eta)*Rd*(3*(Tp-1)+6*(Tp-1)^2*Theta0(eta)+3*(Tp-1)^3*Theta0(eta)^2)/(3*a3*Pr)+4*Rd*(diff(Theta0(eta), eta))^2*(6*(Tp-1)^2*Theta1(eta)+6*(Tp-1)^3*Theta0(eta)*Theta1(eta))/(3*a3*Pr)+4*Rd*(diff(Theta1(eta), eta))*(diff(Theta0(eta), eta))*(6*(Tp-1)+6*(Tp-1)^3*Theta0(eta)^2+12*(Tp-1)^2*Theta0(eta))/(3*a3*Pr)]:

>	(LB,UB):= (0,1): BCs:= [ U0(0) = 0, U1(0) = 0, N0(0) = 0, N1(0) = 0, Theta0(0) = 0, Theta1(0) = 0, U0(1) = 0, U1(1) = 0, N0(1) = 0, N1(1) = 0, Theta0(1) = 1, Theta1(1) = 0 ]:

>

>	Params:= Record( M1= 1.2, Rd=0.8,la=0.8,n1=1.2,Q=0.2,Pj=0.001,Ra=0.8,Ec=1, Pr= 21, delta= 0.2, t1= (1/4)*Pi, lambda0=2,lambda1=3, Qc= 0.1, l1= 1,K=0.4,H=3 ,deltat=0.05 ):

>	NBVs:= [ a1*D(U0)(0) = `C__f` , # Skin friction coefficient (a4+(4Rd(1/3))(1+(Tp-1)(Theta0(0)+0.1e-2exp(l1t1)Theta1(0)))^3)((D(Theta0))(0)+0.1e-2exp(l1t1)(D(Theta1))(0)) = `Nu` # Nusselt number ]: Nu:= `Nu`: Cf:= `C__f`:

>

Solve:= module()
local
   nbvs_rhs:= rhs~(:-NBVs), #just the names
   Sol, #numeric dsolve BVP solution of any 'output' form
   ModuleApply:= subs(
      _Sys= {:-ODEs[], :-BCs[], :-NBVs[]},
      proc({
          M1::realcons:= Params:-M1,
         Pr::realcons:= Params:-Pr,
         Rd::realcons:= Params:-Rd,
         la::realcons:= Params:-la,
         Tp::realcons:= Params:-Tp,
         n1::realcons:= Params:-n1,
         Q::realcons:= Params:-Q,
         Pj::realcons:= Params:-Pj,
         Ra::realcons:= Params:-Ra,
         Ec::realcons:= Params:-Ec,
         t1::realcons:= Params:-t1,
         delta::realcons:= Params:-delta,
         lambda0::realcons:= Params:-lambda0,
         lambda1::realcons:= Params:-lambda1,
         Qc::realcons:= Params:-Qc,
         K::realcons:= Params:-K,
         l1::realcons:= Params:-l1,
         H::realcons:= Params:-H
      })
         Sol:= dsolve(_Sys, _rest, numeric);
         AccumData(Sol, {_options});
         Sol
      end proc
   ),
   AccumData:= proc(
      Sol::{Matrix, procedure, list({name, function}= procedure)},
      params::set(name= realcons)
   )
   local n, nbvs;
      if Sol::Matrix then
         nbvs:= seq(n = Sol[2,1][1,Pos(n)], n= nbvs_rhs)
      else
         nbvs:= (nbvs_rhs =~ eval(nbvs_rhs, Sol(:-LB)))[]
      fi;
      SavedData[params]:= Record[packed](params[], nbvs)
   end proc,
   ModuleLoad:= eval(Init);
export
   SavedData, #table of Records
   Pos, #Matrix column indices of nbvs
   Init:= proc()
      Pos:= proc(n::name) option remember; local p; member(n, Sol[1,1], 'p'); p end proc;
      SavedData:= table();
      return
   end proc ;
   ModuleLoad()
end module:

>

#procedure that generates 3-D plots (dropped-shadow contour + surface) of an expression

ParamPlot3d:= proc(
   Z::{procedure, `module`}, #procedure that extracts z-value from Solve's dsolve solution
   X::name= range(realcons), #x-axis-parameter range
   Y::name= range(realcons), #y-axis-parameter range
   FP::list(name= realcons), #fixed values of other parameters
   {
      #fraction of empty space above and below plot (larger "below"
      #value improves view of dropped-shadow contourplot):
      zmargin::[realcons,realcons]:= [.05,0.15],
      eta::realcons:= :-LB, #independent variable value
      dsolveopts::list({name, name= anything}):= [],
      contouropts::list({name, name= anything}):= [],
      surfaceopts::list({name, name= anything}):=[]
   }
)
local
   LX:= lhs(X), RX:= rhs(X), LY:= lhs(Y), RY:= rhs(Y),
   Zremember:= proc(x,y)
   option remember; #Used because 'grid' should be the same for both plots.
      Z(
         Solve(
            LX= x, LY= y, FP[],
            #Default dsolve options can be changed by setting 'dsolveopts':
            'abserr'= 0.5e-7, 'interpolant'= false, 'output'= Array([eta]),
            dsolveopts[]
         )
      )
   end proc,
   plotspec:= (Zremember, RX, RY),
   C:= plots:-contourplot(
      plotspec,
      #These default plot options can be changed by setting 'contouropts':
      'grid'= [25,25], 'contours'= 5, 'filled',
      'coloring'= ['yellow', 'orange'], 'color'= 'green',
      contouropts[]
   ),
   P:= plot3d(
      plotspec,
      #These default plot options can be changed by setting 'surfaceopts':
      'grid'= [25,25], 'style'= 'surfacecontour', 'contours'= 6,
      surfaceopts[]
   ),
   U, L #z-axis endpoints after margin adjustment
;
   #Stretch z-axis to include margins:
   (U,L):= ((Um,Lm,M,m)-> (M*(Lm-1)+m*Um, M*Lm+m*(Um-1)) /~ (Um+Lm-1))(
      zmargin[],
      (max,min)(op(3, indets(P, 'specfunc'('GRID'))[])) #actual z-axis range
   );
   plots:-display(
      [
         plots:-spacecurve(
            {
               [[lhs(RX),rhs(RY),U],[rhs(RX),rhs(RY),U],[rhs(RX),rhs(RY),L]], #yz backwall
               [[rhs(RX),rhs(RY),U],[rhs(RX),lhs(RY),U],[rhs(RX),lhs(RY),L]] #xz backwall
            },
            'color'= 'grey', 'thickness'= 0
         ),
         plottools:-transform((x,y)-> [x,y,L])(C), #dropped-shadow contours
         P
      ],
      #These default plot options can be changed simply by putting the option in the
      #ParamPlot3d call:
      'view'= ['DEFAULT', 'DEFAULT', L..U], 'orientation'= [-135, 75], 'axes'= 'frame',
      'labels'= [lhs(X), lhs(Y), Z], 'labelfont'= ['TIMES', 'BOLDOBLIQUE', 16],
      'caption'= nprintf(cat("%a = %4.2f, "$nops(FP)-1, "%a = %4.2f"), (lhs,rhs)~(FP)[]),
      'captionfont'= ['TIMES', 14],
      'projection'= 2/3,
      _rest
   )
end proc:

>

>	ParamPlot3d( GetNu,Q= 0..5, Rd= 0..5, [ Pr= 21 ], labels= [Q, gamma, Nu] );

Error, (in plot/iplot2d/levelcurve) could not evaluate expression

Download P6_3D_plots.mw

How do I plot binned data in Maple?...

Asked by: Traruh Synred 55

September 15 2023

0 6

I have data that I've binned in list. I'd like to plot in as a histogram.

dataplot sort of does it but doesn't give the x-axis that I used but just the bins.

Yes, I know about Histogram from Statistics. To plot 1000000 values I have to enter them all into a list which seems crazy when all I want is a 100 bin histogram. Doing the binning is trivial, but I can't figure out how to plot it with a sensible x-axis (show range used to define the 'histogram')

Nothing works after updating from Maple 2015 to 20...

Asked by: leouam 25

September 15 2023

0 3

Hello,

I am experiencing difficulties using my old Maple programs with the newer version. I tried changing the types of inputs and the typesetting level, but it just doesn't work. I would appreciate it if someone could help me overcome my ignorance.

Some simple input is attached with the output.

how to input complex function in a complex differe...

Asked by: bashar27 55

August 20 2023

0 0

restart;
alias(u = u(x, y, t), f = f(x, y, t));
u, f
u := (f+sqrt(R))*exp(I*R*t);
/ (1/2)\
\f + R / exp(I R t)
pde1 := I*(diff(u, t))+diff(u, x, x)+2*lambda*u*abs(u)*abs(u)-gamma*(diff(u, x, t));
// d \ / (1/2)\ \
I ||--- f| exp(I R t) + I \f + R / R exp(I R t)|
\\ dt / /

/ d / d \\
+ |--- |--- f|| exp(I R t)
\ dx \ dx //

/ (1/2)\ 2
+ 2 lambda \f + R / exp(I R t) (exp(-Im(R t)))

2
| (1/2)|
|f + R |

// d / d \\ / d \ \
- gamma ||--- |--- f|| exp(I R t) + I |--- f| R exp(I R t)|
\\ dx \ dt // \ dx / /

DTM method with Dsolve method comparision...

Asked by: SHIVAS 35

August 18 2023

0 0

how to plot graphs for both methods and comparison of different method values for Diff(f(eta),eta, eta) at eta =0

>

NULL

>

for k from 0 to n[1] do F[k+4] := solve((1+a)*(k+1)*(k+2)*(k+3)*(k+4)*F[k+4]-a*(k+1)*(k+2)*G[k+2]-R*(sum(F[k-m]*(m+1)*(m+2)*(m+3)*F[m+3], m = 0 .. k))+R*(sum((k-m+1)*F[k-m+1]*(m+1)*(m+2)*F[m+2], m = 0 .. k)), F[k+4]) end do

-(1/60)*(R^2*a2*a3*al-3*R^2*a4*al^2+2*R*a*a3^2-R*a*al*G[2]+2*R*a3^2-3*a^2*G[3]-3*a*G[3])/(1+a)^2

(1)

>

-(1/6)*(R^2*a2*a5*al*c^2-R^2*a6*al^2*c^2+2*R*a*a3*al*c-2*R*a*a5*al*c+2*R*a3*a5*b*c+6*a*a4*b-2*a*a6*b)/b^2

(2)

>

-p3*Q[4]-(1/24)*p1^2*a2^2*a7+(1/24)*a2*p1^2*al*a8-(1/12)*p1*a2*p3*Q[2]-(1/12)*p1*a3*a8-(1/4)*p1*a4*a7-(1/24)*a2*a7*al^2*p1^3+(1/24)*a8*al^3*p1^3-(1/12)*al^2*p1^2*p3*Q[2]-(1/12)*al*a3*a7*p1^2-(1/4)*p1*al*p3*Q[3]

(3)

>

(1/6)*p2*(a10*al^2*p2-a2*a9*al*p2-2*a3*a9*p4-2*a3*a9)/(p4+1)^2

(1/24)*p2*(a10*al^3*p2^2-a2*a9*al^2*p2^2+a10*a2*al*p2*p4-a2^2*a9*p2*p4-2*a3*a9*al*p2*p4+a10*a2*al*p2-2*a10*a3*p4^2-a2^2*a9*p2-2*a3*a9*al*p2-6*a4*a9*p4^2-4*a10*a3*p4-12*a4*a9*p4-2*a10*a3-6*a4*a9)/(p4+1)^3